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SUMMAEY 


.  The  laminar  flow  of  a  slightly  viscous  Incompressible  fluid  that 
issues  from,  a  slit  and  passes  over  a  flat  plate  is  investigated  in  a 
region  far  enou^  from  the  slit  for  the  boundary— layer  equations  to  be 
valid.  By  assuming  similar  velocity  distributions  along  the  plate,  the 
partial,  differential  equation  for  the  boundary  layer  is  reduced  to  a 
third-order  nonlinear  ordinary  differential  equation.  This  equation  is 
integrated  by  nmerical  means  for  the  required  boundary  conditions.  The 
solution  gives  the  velocities  at  points  in  the  fluid  and  the  sinrface 
friction  .at  points  on  the  plate. 

Some  of  the  specific  results  obtained  are  that  the  velocities  paral¬ 
lel  tp  the  plate  vary  inversely  as  the  1/2  power  of  the  distance  from 
the  slit,  that  the  velocities  perpendicular  to  the  plate  vary  inversely 
as  the  3/^  of  the  distance  from  the  slit,  that  the  width  of  the 

disturbed  region  Increases  as  the  3/^  power  of  the  distance  from  the 
slit,  that  the  surface  friction  is  independent  of  the  viscosity  of  the 
fluid  and  varies  inversely  as  the  5/^  power  of  the  distance  from  the 
slit,  that  the  rate  at  which  moment imi  parallel  to  the  plate  passes 
throu^  a  plane  normal  to  the  plate  varies  Inversely  as  the  l/h  power 
of  the  distance  from  the  slit  to  the  normal  plane,  and  that  the  quantity 
of  fluid  passing  throu^  a  plane  normal  to  the  plate  per  unit  time  varies 
directly  as  the  l/h  power  of  the  distance  from  the  slit  to  the  normal  plane 


INTEODUCTION 


The  laminar  flow  of  an  infinite  uniform  stream  of  incompressible 
fluid  of  small  viscosity  over  a  flat  plate  at  zero  angle  of  attack  was 
first  investigated  analytically  by  Prandtl  (reference  1).  By  a  proper 
choice  of  variables  the  partial  differential  equation  of  the  Prandtl 
boundary  layer  was  reduced  to  a  norxlinear  ordinary  differential  equation 
of  the)  third  order.  An  accurate  solution  of  the  equation  that  gave  the 
skin^friction  coefficient  for  the  plate  and  the  velocities  in  the  field 
of  flow. was  first  obtained  by  Blasius  (reference  2)  and  later  by  others. 
(For  example,  see  references  3  to  6.)  The  nimierlcal  values  obtained  by 
the  different  investigators,  all  of  whom  used  either  series  or  numerical 


methods  to  integrate  the  nonlinear  differential  equation,  were  about  the 
same . 

The  two-dimensional  laminar  flow  of  a  slightly  viscous  incompressible 
fluid  from  a  slit  into  an  infinite  region  of  still  fluid  was .investigated 
analytically  by  Schlichting  (reference  7)  and  by  Bickley  (reference  8). 

By  introducing  new  variables  into  the  Prandtl  boundary-layer  equation 
and  by  using  the  fact  that  the  rate  of  flow  of  momentum  through  a  cross 
section  of  the  Jet  is  independent  of  the  distance  from  the  slit,  the 
partial  differential  equation  of  the  boundary  layer  was  reduced  to  a 
third— order  nonlinear  ordinary  differential  equation.  Schlichting  inte¬ 
grated  the  equation  by  the  method  of  series  and  Bickley  later  Integrated 
the  equation  in  a  closed  form.  The  solutions  gave  the  velocities  at 
points  in  the  field  of  flow.  In  addition  to  the  solutions  for  the  flow 
over  a  fiat  plate  and  for  the  flow  from  a  slit,  reference  9  also  gives 
the  few  other  known  solutions  of  the  boundary— layer  equations.  The  solu¬ 
tions  are  given  for  the  flows  near  the  stagnation  points  of  plane  bodies  • 
and  of  bodies  of  revolution,  for  the  flow  along  a  wall  in  a  converging 
channel,  and  for  the  flow  in  a  round  Jet  from  which  the  fluid  issues 
from  a  small  hole  in  a  wall.  Any  additional  solution  of  the  boundary- 
layer  equations  is  then,  aside  from  any  practical  application,  of  some 
Importance  in  boundary— layer  theory. 

The  purpose  of  the  present  work  is  to  investigate  by  means  of  the 
Prandtl  boundary— layer  equations  the  laminar  flow  of  a  slightly  viscous 
incompressible  fluid  that  issues  from. a  slit  and  passes  over  a  flat 
plate.  The  arrangement  (fig.  l)  may  also  be  interpreted  as  the  laminar 
flow  into  a  very  large  container  of  still  fluid  from  a  slit  at  the 
intersection  of  two  of  the  walls.  By  use  of  the  Prandtl  boundary-layer 
equations  together  with  the  momentum  theorem,  a  substitution  involving 
other  variables  was  found  for  the  variables  that  appear  in  the  Prandtl 
boundary— layer  equation.  The  substitutions  reduced  the  boundary— layer 
equation  from  a  partial  differential  equation  to  a  nonlinear  ordinary 
differential  equation  of  the  third  order.  The  equation  was  Integrated 
numerically  for  the  req'uired  boundary  conditions  with  the  aid  of  the 
general  purpose  computing  system  of  the  Bell  Telephone  Laboratories  in 
the  Langley  Bell  computing  section. 

The  reduction  of  the  boundary— layer  partial  differential  equation 
to  an  ordinary  differential  equation  in  the  cases  of  the  flat  plate,  the 
Jet,  and  the  combined  flat  platd  and  Jet  flow  is  made  by  finding  suitable 
substitutions  for  the  variables  that  appear  in  the  Prandtl  boundary-layer 
equations.  By  making  these  substitutions  the  velocity  distributions  in 
planes  nonnal  to  the  direction  of  the  main  flow  are  assumed  to  be  similar 
The  success  of  the  substitutions  in  reducing  the  partial  differential 
equation  of  the  boundary  layer  to  an  ordinary  differential  equation  in 
these  special  cases  means  that  the  assumption  of  similar  profiles  is 
compatible  with  the  equations  describing  the  motion  and  that  the  similar 
profiles  exist  wherever  the  assumptions  made  in  deriving  the  equations 
that  describe  the  motion  are  valid. 
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The  fact  that  the  profiles  are  the  similar  profiles  predicted  by 
the  laminar-boundary-layer  theory  has  been  verified  experimentally  both 
for  the  flow  of  a  uniform  Infinite  stream  over  a  flat  plat©  (reference  10) 
and  for  the  flow  In  a  round  jet  (reference  11).  The  experimental  Investi¬ 
gation  of  the  flow  from  the  Jet  verified  the  analysis  (reference  9)  for 
all  regions  except  In  the  Immediate  neighborhood  of  the  orifice.  The 
application  of  the  Prandtl  boundary— layer  eq^uatlons  to  the  present  case, 
an  application  that  Implies  assumptions  similar  to  those  made  in  analyzing 
the  flow  of  a  uniform  stream  over  a  flat  plate  and  in  analyzing  the  flow 
from  a  round  orifice,  is  therefore  expected  to  lead  to  results  that  are 
valid  for  regions  not  too  close  to  the  slit. 


SYMBOLS 


y 


nondime  ns  ional  distance  parallel  to  siarface  of  plate  (5c/a) 
distance  parallel  to  surface  of  plate 

nondimens ional  distance  normal  to  surface  of  plate  (y/a) 
distance  normal  to  surface  of  plate 
nondimensional  velocity  component  in  x-direction 
velocity  component  in  x-4irection 


nondimensional 


velocity  component  in  y— direction  (v/^a) 


velocity  component  in  y— direction 
density 

coefficient  of  viscosity 
kinematic  viscosity 

reference  velocity 
reference  length 

reference  Reynolds  nimiber  (u^a/v^ 
function  of  j] 


function  of  x 


f -unction  of  x 


K  =  —  = 

°Sa  r  2 

f®d,| 


static  pressure 


constant 


Cj,  constants  of  Integration 


momentum  in  x— direction 


reference  momentum  in  x— direction,  one  side  of  plate 


lower  limit  of  |  in 


G  =  l 


a  = 


surface  friction  f  \i 


quantity  of  flow,  one  side  of  plate 


maximum  velocity  at  a  section 

^max  distance  from  saarface  of  plate  to  point  of  maximum  velocity 
Sutscripts: 

0  conditions  at  surface  of  plate 

00  conditions  at  infinity 


ANALYSIS 


After  assuming  that 
pressure  gradient  exists 
the  usual  assumptions  of 
the  eq^uations  describing 


motion  vith 


^  =  0 
5x 


the  flow  is  incompressible j  that  no  static- 
in  the  direction  of  x-axis  (fig.  l)^  and  that 
the  boundary— layer  theory  (reference  9) 
the  motion  are:  The  boundary— layer  eq^uation  of 


u 


Sx 


+  V 


(1) 


the  eq^uation  of  continuity  for  incompressible  flow 


^+^=0  (2) 
Sx  hy 


and  the  momentum  theorem  applied  to  the  x— component  of  the  flow 


noo 


(Jo 


pn^dy 


dx 


(3) 


In  order  to  make  the  eq^uations  nondimensional,  the  velocities  are  divided 
by  a  reference  velocity  and  the  lengths  by  a  reference  length.  The 
nondimensional  velocities  and  lengths  are 


Now  eq^uation  (3)  becomes 


V  Z- 


X 

X  =  — 


d  dy 


^a  (  hu 


d  u  dy 

Jo 


1  f'du. 


^a\Sy/o 


where 


2^0  =  — ^ 
a  y 


.  The  velocity  distributions  in  planes  perpendicular  to  the  x-^xis 
are  now  assumed  to  be  similar  to  one  another.  Velocity  distributions 
are  said  to  be  similar  if  all  the  velocity  distributions  are  given  by 
one  curve  when  u/u^^y  is  plotted  against  y /y^^y .  In  order  to  make 

the  assumption  of  similar  profiles,  the  independent  variables  y  and 
are  replaced  by  new  independent  variables  t)  and  x  by  means  of  the 
following  substitutions  (reference  12); 


•  ••• 
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where  f  =  f(T|)  and  f  =  0(x),  and 

•  •  • 

•  •  • 

•  •• 


where  P  =  P(x).  The  expressions  for 
fixed  value  of  x  then  are 


u/u 
/  max 


and  for 


y/y. 


max 


for  a 


and 


% 


u 

%iax 


_y _  ^  “H 


All  the  velocity  distributions  in  planes  perpendicular  to  the  x-axis  are 
therefore  given  by  one  curve  of  u/u^^  plotted  against  y/y^nax 

consequently  the  foregoing  substitutions  for  y-  and  x  mean  that  similar 
profiles  have  been  assumed.  Whether  the  assumption  of  similar  profiles 
is  compatible  with  equations  (l)  to  (3)  which  describe  the  motion  is 
decided  by  whether  these  equations  reduce  to  ordinary  differential 
equations. 


Determination  of  Relations  for  0,  and  u 

The  functions  0  and  p  are  to  bo  determined  by  using  equations  (l) 
to  (3).  If  the  values 


p 


dT) 


and 


5u 


•  ••• 
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are  substituted  in  equation  the  resiiLt  is 


•  •  • 

•  • 
•  • 

•  '  • 


0^  u 


Ooo 


0 


2 

f  dq 


dbc 


1_  1 

01 


Uti/q 


or,  since 


n  00 


f  cLt]  is  independent  of  x 


Nov  let 


then 


i  M  _  i  M 

P  cbc  0  dx 


^Tl/O 


L 


noo 


2  P‘ 

f  dq 


liO 


P  dx  0  dx  p2 


(5) 


Equation  (5)  can  be  solved  by  using  equation  (l)  to  provide  a  relation, 
which  does  not  Involve  x,  between  0  and  p.  After  equation  (l)  is 
made  nondimensional,  the  result  is 


# 


^  Su  ^  Su  _  1 
Sx  5y  5^2 


(6) 
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Eq,uation  (6)  is  to  te  rewritten  in  terms  of  f,  P,  and  t) 

the  relation 


•  ••• 
•  ••• 


•  •  • 
•  •• 
•  • 


f 

u  =  — 
0 


it  follows  that 


and 


Su  _  ^  df 
Sy  0P  dT) 


1  d^f 


0p^  dT]^ 


% 


/  M. 

5u  _  I  dtj  dbc  f  d0 
^  \  ^ 


From  eq^uatlon  (2) 


V  = 


0 


'bsr 


py 


dy  = 


JO 


5u 

5x 


dy  =  - 


Jo 


3  dT] 


where 


or 


dy  =  3  dT) 


V  = 


«  pr,  p  M 

ta  df  ,  db: 

—  r\  —  dt]  +  — — 

0  Jo  0^ 


f  dTi 


.  From 
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By  partial  integration 


•  ••• 

•  ••• 

•  •  • 
•  •  • 
•  •• 


df 


T]  —  dT]  =  fr|  —  f  dr) 

0  Jo 


then,  the  result  for  v  is 


0  Jo  ^02  dx  0  dx 


If  values  for  u,  v,  — ,  — ,  and  are  substituted  in  eq.Tia- 

Bx^  ay’  Sy2 

tion  (6),  the  result  is 


df 


f  (  ^T)  f_ 

0  \  P0  dx  02  dx^ 


jiM.r  f 

0  Jo  ^02  0  dx^ 


=  (7) 

0p  dTi  Ea  0p2  dTi2 


The  left-hand  and  rl^t— hand  sides  of  equation  (7)  are  multiplied 
hy  0^,  and  the  factor  0/p^  which  then  appears  on  the  right-hand  side 
of  the  equation  is  replaced  hy  the  value  of  0/p^  from  equation  (5)* 
After  further  simplification,  equation  (7)  can  he  written  as 


1  ^ 
0  dx 
1  ^ 
3  dx 


dq 


nq 


Jo 


f  dq  -  —  — 
I^aK  dT)2 


.f2  +  ^ 
d7 


f  d  2  d^f 

''  RaK  dti^ 


(8) 


The  variables  x  and  q  are  independent.  Consequently,  a  change  in  x 
does  not  change,  the  right-hand  side  of  equation  (8)  and  therefore  does 
not  change  the  left-^iand  side  either.  A  change  in  q  likewise  does  not 
H  *  change  the  left^and  side  of  equation  (8)  and  thus  does  not  change  the 

right-hand  side  either.  Both  the  left-hand  and  right-hand  sides  of  equa¬ 
tion  (8)  are  therefore  equal  to  a  constant  C.  Then 

« 


f 
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•  ••• 
•  ••• 

•  •  • 
•  •  • 
•  •• 
•  • 


or 


or,  by  integrating 


1  djZi 

0  dx 

1  dP 
(B  dx 


=  C 


<10  _  dp 

r  r 


i*  -  CiP® 

Vhen  eq[uation  (9)  is  substituted  in  equation  (5),  the  result  is 


pi-c  ^  ^ 


dx  1  -  2C 


or,  after  integrating  ^or  ^  ^  2^ 


P  = 


KCi 


-ft-  f 

2C  -  1  2 


X  +  C_  (2  -  C) 


1 

2-C 


i 


and 


Then 


u  = 


0  =  C; 


KC 


2C  -  1 


X  +  (2  -  C) 


C 

2-C 


f  ^ 


/  \ 

c 

# 

Cl 

2-C 

1 

KC 


- 'X  +  Cp  I  (2 

2C  -  1 


-c) 


2-C 


(9) 
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The  expression  for  u  contains  two  arbitrary  constants  and  Cg. 

The  constant  is  a  scale  factor;  the  constant  C2  is  a  translation, 

factor  that  determines  the  location  of  the  slit  along  the  x-axis.  There¬ 
fore,  the  constant  can  be  made  unity  and  the  constant  C2  can  bo' 

made  zero  without  affecting  the  generality  of  the  solution  in  any  way 
other  than  fixing  the  origin  of  the  slit  at  x  =  0.  The  expressions 
for  P,  0,  and  u  therefore  become 


P  = 


Kx 


2C  -  1 


(2  -  C) 


1 

2-C 


0  = 


Kx 


2C  -  1 


(2  -  C) 


C 

2-C 


and 


u 


Kx 


2C  -  1 


(2  -  C) 


C 

2-C 


1. 


Kx 


2C  -  1 


(2  -  C) 


1 

2-C 


The  constant  C  cein  he  evaliiated  hy  use  of  the  expression  for  the 
momentum  in  the  x— direction: 


M  =  p  I  u^dy 


or 


poo 


M  /  2,  P 

_  .  u  dy  =  — 

pUa^a  Jo  f 


f^dT| 


or 


•  ••• 
•••• 


•  •• 
•  • 


i'3 


M 

pUa^a 


Kbc 

2C  -  1 


1- 2C 

2- C 


1  00 

f^dT) 

Jo 


(10) 


The  arbitrary  length  a  is  interpreted  as  the  distance  along  the 
x-axis  from  the  origin  to  the  point  where  the  similar  profile  has  been 
established.  At  this  distance  a  from  the  origin  the  rate  of  momentum 

flow  is  equated  to  pUg^  a;  thus,  the  reference  velocity  Ug  is  defined 
When  x=a  orx=l,  M=Mg  and,  consequently,  equation  (10)  at 
X  =  a  becomes 


% 


1- 2C 

2- 0 


poo 

r^dT) 


(Jo 


(11) 


The  expression 


0  00 

f^dq 

Jo 


is  a  constant. 


Let 


f^dT)  =  1 

Jo  ' 


The  fmetion  f  to  be  found  must  satisfy  this  condition.  Then  equa- 
tion  (11)  hecomes 


^  •  Therefore 


1- >2C 

2- C 


2K  +  1 


0 


C 


K  +  2 


•  ••• 


t 


Ik 


Then 


•• 


••  • 
•  •  • 
•  •• 


(12) 


and 


2K+1 

0  =  X  ^ 


(13) 


and 


X  >  1 


It  can  now  he  shown  that  the  requirement  that  C  ^  2  and  that 


which  was  stipulated  in  determining  P  as  a  function  of  x  by  using 
equations  (5)  and  (9)^  Is  satisfied.  The  expression  that  was  previously 
obtained  for  C  is 


C 


2K  +  1 
K  +  2 


From  physical  considerations  0  <  (  — ^  <  oo  and  0  <  <  oo.  There— 

fore,  0  <  K  <  Thus,  because  K  >  0,  the  expression  for  C  gives 


C  >i 

2 


4 


•  ••• 
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Because  K  < 


C  <  2 


•  • 


Therefore, 


C  /  2 


and 


In  order  to  shov  that  the  flow  appears  to  come  from  a  slit,  the 
edge  of  the  region  in  which  the  velocity  u  is  other  than  zero  can  he 
arbitrarily  defined  hy  choosing  a  fixed  value  of  t)  that  corresponds  to 
a  small  fixed  value  of  u.  Thus,  let  the  edge  of  the  jet  region  he 
given  hy 


4 


T) 


^edge 


=  Constant 


Then 


'edge 


^edge 

^+2 

3 


=  Constant 


or 


K+2 


Jedge  =  ^  ^  X  Constant 


% 


The  flow  therefore  does  appear  to  come  from  a  slit  located  at 

(x  =  0,  y  =  0),  The  slit  can  he  placed  at  a  position  of  x  that  is 

not  zero  hy  choosing  a  value  for  Cg  that  is  not  zero.  When  the  flat 

\ 

plate  is  removed,  (  —  ]  =0  or  K  =  0;  and  then  from  equation  (Ik) 

\^Tl/0 


l6 


u  =  f( 

1/3  I  2/3 


which  is  similar  to  the  relation  for  u  in  reference  8. 


Derivation  of  the  Differential  Eq,uation  for  G 

The  first  step  in  determining  the  differential  eq.uation  for  G  is 
to  suhstltute  equations  (12)  and  (I3)  in  equation  (7).  The  equation  of 
motion,  therefore,  hecomes 


2K  +  1^2  ^ 


L-  ^1. 

^a  dq^ 


The  assumption  of  similar  profiles  has  reduced  the  partial  differential 
equation  of  the  hoimdary  layer  to  an  ordinary  differential  equation; 
thus,  an  indication  is  given  that  the  assumption  of  similar  profiles  is 
consistent  vith  equations  (l)  to  (3)  which  describe  the  motion.  Now, 


f  dT]  =  F 


where  F  =  F(q).  Then,  because 


d^f  ^  ^ 
dq^  dq^ 


equation  (I7)  iDecoiries 


•  ••• 

•••• 

•  •  • 
•  •  • 
•  •• 


1  ^  1  -  K  j,  d^F  SK  +  1 

d,3  3  4,"=  *  3 


Nov,  let 


gq 


and 


7  = 


¥  =  yjr 


where  G-  =  G-(|)«  Thus, 


and 


dE  gy  dG 

dT]  ^  3 


d^g  ^  g-^y  d-% 
dT|^  9  d,|2 


Equation  (l8)  then  'becomes 


^3p  ^2- 

+  (1  _  K)G  ^  +  (1  +  2K) 
d|3  d|^ 


8 


where 


•  ••• 
•  •  •  • 


#  •  • 
•  •• 
•  • 


K  = 


dfN 


d%' 


R 


a 


\d|‘ 


0 


9Ea^  Vi^Jo 


Nov,  let 


a  =  9^/^R„^/^ 


Then,  equation  (I9)  becoiries 


d^G 

d|3 


d^G 

^e2 

.d|  /Q 


dl^ 


1  +  2 


d% 


d|' 


'0 


.  0 

d|j 


(20) 


The  conditions  to  he  satisfied  hy  G  are,  at  ^  =  0, 


G  =  0 


dll  " 


and,  at  |  =  00, 


G  =  Constant 


n 


Integration  of  the  Differential  Equation  for  G 


The  differential  equation  for  G,  equation  (20),  is  integrated 
numerically  hy  use  of  the  method  of  reference  I3.  The  numerical  inte¬ 


gration  vas  started  hy  computing  the  values  of  G, 


and 


for  the  first  five  intervals  of  ^  hy  means  of  a  Taylor^  s  series.. 
Taylor^s  series  that 'is  valid  for  small  values  of  i  is 


The 
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From  the  nature  of  the  problem  each  term  of  the  eq-uatjon  Is  continuous 
in  the  interval  from  0  to  oo;  the  equation  may  therefore  he  integrated 
directly  term  hy  term.  Then 


••  f 
•  •  • 
m  •• 
•  • 


noo 


G  ^  d|  + 
dl^ 


P 

_  W7o_ 

Jc 

G^^dl. 


1  +  2 


d2G 

/Q 


r)oo 


0 


<f) 


Integrating 


d3G 

Gr  -  d|  hy  parts  gives 

4|^ 


Jo 


d|3 


fG 

d|2 


0  liO 


d^G  dG 
—  —  d| 

d|2  d| 


n“  -2 


do 


d^G  ^ 
d|2  d| 


Integrating 


poo  ^ 

d^g  ^ 

Jo  4|24| 


d|  hy  parts  gives 


d%  dO 
to  ai=  '15 


■'S 


00  f>  00 


dl2 


dl  =  - 


d^G  dG 


(O  d| 


2  dl 


dl 


Therefore 


00  P 

O^d|  =  o 


ncc  3 


G  —  dl  =  0 

Jo 


and  thus 
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Nov,  irxtegrating 


^2  jfc  . 

G-  — ^  loj  parts  gives 

dr 


•  •  • 
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pOC  00 

£s  di  .  o2 

d-5 

0  '^5  Jn 


^  2G  ^  d| 
d|  d| 


Then  equation  (21)  'becomes 


-2  1 


+  1  +  2 


d|  =  0 


The  quantity  G  is  always  positive,  except  on  the  plate  where  it 
is  zeroj  therefore,  equation  (22)  cannot  be  time  unless  -  i. 


Consequently, 


is  taken  as  the  starting  value  of 


the  numerical  integration  of  equation  (20). 

The  equation  was  integrated  hy  use  of  the  ’’fourth  approximation”  of 
reference  13  for  an  increment  in  r  =  0.10;  the  numerical  integration 
"began  at.  ^  =  O.5O.  In  order  to  determine  whether  the  inoreiaent  in 
I,  A|  =  0.10,  was  sufficiently  small,  the  equation  was  also  integrated 
for  an  increment  in  A|  =  0.05;  the  numerical  integration  began 
at  I  =  0,50.  Some  of  the  results  from  the  integration  with  A|  =  O.O5, 
the  more  accurate  integration,  are  given  in  table  I.  A  comparison  of 
some  of  the  values  of  G  and  its  derivatives  o'btained  by  integrating 
with  A|  =  0.05  and  with  A^  =  0.10  is  given  in  table  II;  the  values 
at  ■  I  =  0.60  computed  by  the  Taylor’s  series  are  also  included  for 
comparison  with  those  computed  by  Falhner’s  method  (reference  I3).. 
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The  results  of  the  two  integrations  agree  to  five  significant 

figures  up  to  |  =  8,50  at  which  value  —  is  less  than  1  percent  of 

its  maximum  value.  From  |  =  8,50  to  |  =  I3.IO,  at  which  value  the 
computation  for  A|  =  0.10  was  stopped,  the  agreement  between  the 
computations  for  the  two  different  increments  in  |  becomes  poorer. 

The  values  of  0,  however,  are  the  same  to  five  significant  figures  for 
both  computations.  At  ^  =  I3.OO,  0  has  already  attained  its  asymptotic 
value,  a  value  that  is  necessary  for  the  computation  of  v  at  |  =  00, 

(See  table  I.) 

As  a  further  check  on  the  accuracy  of  the  nvimerical  integration, 
the  following  condition  was  used: 


T  00 

2 

f  dri  =  1 

Jo 

which  is  the  same  as 

4 


(23) 


The  values  of 


PI 


(JO 


obtained  by  using  Weddle’s  rule  (reference  Ik) 


vith  A|  =  0,05  are  given  in  table  III.  It  is  seen  that  the  foregoing 
condition  (equation  (23))  is  confirmed  to  five  significant  figui’es.  In 


order  to  determine  whether  the  integral 


vith  sufficient  accuracy  by  the  integral 


noo  Q 

0 


112.30 

/ 

^dG\2 

( 

Jo 

d|  is  approximated 


)  d|,  the  assumption 


vas  made  that  for  large  values  of  |  the  curve  of 
could  be  approximated  by  the  function 


plotted  against 


I 


^  =  Ae-®^ 


•  ••• 
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The  constants  A  and  B  were  evaluated  In  the  range  of  |  from  |  =  9-^5 
to  I  =  114-.95*  The  constant  B  was  observed  to  Increase  slowly.  For 


•  •  • 
•  •  • 
•  •• 
•  • 
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dl 


the  Integral  from  to  oo  is 


noo  p 

3£fd|  = 


2B 


(0.00009149710)^ 

2(0.9452) 


=  0.4425  X  10 


-8 


where 


=  12.30 


^  =  0.00009145710 


and 


B  =  0.9452 


Because  —  decreases  as  It  increases  and  because  B  Increases 
d| 

fi  “  _ 2 


as  Ij^  Increases,  the  Integral 


dG  \ 

—  \  d|  is  given  accurately  to 


five  significant  figures  by  the  integral 


pis. 30 


Jo 


<iiy 


In.  order  to  determine  whether  the  integral 


mated  with  sufficient  accuracy  hy  the  integral 


0 

nh.95 


d I  is  appro:-:!— 


do 


d| 


d|,  the  same 


dG 

fimctlon  for  —  for  large  values  of  f  was  assumed:  that  is,  for 
d|  ^  ’  ’ 


if 


the  integral  from  to  oo  ig 


(f 


n“  \d|A 

^  ^ 

d|  B 


where 


and 


=  0.000007^33791 
0.950682 

=  0.00000782 


d| 


lh.95 

0.000007^33791 


i 


B  =  0.950682 
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dG 

Because  —  decreases  as  increases  and  Because  B  increases 

d? 

noo 


as  increases,  adding 

l/11^.95 

does  not  change  its  value  to  five  significant  figures. 


—  d|  to  the  value  of  G  at  |  =  lU.95 


noo 


The  last  value  in  table  III  for 


Jo 


d|  and  the  asymptotic 


value  of  G  are  therefore  both  unchanged  by  extension  of  the  curves 
to  I  =  00,  The  results  of  the  computation  for  A|  =  O.O5  are  conse¬ 
quently^  taken  as  correct  to  five  significant  figures  for  0  ^  |  ^8.50 

and  the  asymptotic  value  of  G  is  taken  as  1.2599.  The  curve  of  — 
plotted  against  ?  is  shoTm  in  figure  2. 


Determination  of  Final  Expressions 

Now  that  the  values  of  G,  have  been  tabulated, 


d^3 


the  completion  of^the  solution  of  the  problem  requires  that  expressions 
for  u,  u,  V,  V,  Q^,  M,  and  ^  be  determined. 

In  order  to  obtain  the  expressions  for  u  and  u,  use 


f 

u  =  — 

0 


2K+1 

3 


3  dg  ^  9^/3p  1/3  Ji 


xl/2 


xl/2 


X 


(24) 


The  expression  for  Is 


E.  = 


_  a 
V  “*  V 


% 

pa 


pv^ 


therefore,  u  may  now  be  written  as 


u 


^  92/3/Maa\V^  f 


3  Vpv^ 


xL/2 


(25) 
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•  The  expression  for  u  is 

•  • 


or 


u 


9^/3  ^  ^ 

3  [p^) 


(26) 


The  u  component  of  the  velocity  therefore  varies  inversely  as  the 
1/2  power  of  the  distance  from  the  slit. 

The  expression  for  v  is 


where 


V  = 


£ll  IP 

0  dx 


f  dr) 


I/O 


L.  M  _  i  iP' 

2  dx  0  dx 


0 


^  CLy  dG 
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p  =  X  3 
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f  dTi  =  7G 
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and 


Then 


at] 

3 


V  =  -L_  IL±^  I  ^  +  2-^  g] 

£t2  V  3  d|  3  j 

X  3 

and,  after  the  appropriate  substitutions  for  y  and  K  have  been  made. 


V  =  9^^^  1  /3  f  g' 

1/3  U  d|  4 


or 


V  = 


9V3  1 
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/MaaV 

i  2  ) 

\PV  / 


-\l/6  ^3/4  V-^  d| 


3?  ^-G 


(27) 


The  expression  for  v  is 


T  =  vUa  = 


9I/3 
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Maa 


pv 


_  e  dG  p 

a/6  _3/4  d| 
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.  9l/3/M,v\l/3  JA2  .  ^  . 
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(28) 
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The  magnitude  of  the  vertical  velocity  v  thus  varies  inversely  as  the 
3/4  pover  of  the  distance  from  the  slit.  Whether  the  vertical  velocity 

is  away  from  or  toward  the  plate  depends  on  the  sign  of  the  term  —  —  G 

in  ec[uation  (28). 


When  I  — >  oOj  it  can  he  shown  that 
as  I — >-00.  Therefore, 


because 


Constdnt 


or 


91/3 


(1.2999)9^3  /MaVV^/3  ^/12 

X 


(29) 


The  flow  from  the  slit  thus  induces  a  velocity  towards  the  plate  at 
sufficiently  large  distances  above  or  below  the  plate. 

The  quantity  of  fluid  passing  through  an  imaginary  plane  normal  to 
the  plate  and  extending  to  infinity  in  one  direction  is  given  as 


Q  = 

^00 


noo 


n  «> 


u  dy  =  U„a 


iiO 


u  dy 


liO 


and,  after  the  appropriate  suhstl  tut  Iona  for  u,  and  y  have  been 

made. 


or,  finally. 


Q, 


0  00 


d| 


•1 
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(30) 
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vhere 


(>00  =  1-2599 


The  quantity  of  fluid  passing  through  a  plane  normal  to  the  plate  per 
unit  time  thus  varies  directly  as  the  1/4  power  of  the  distance  from 
the  slit  to  the  normal  plane. 

The  expression  for  the  surface  friction  t  Ig 


/3u' 


'0 


a  I  Sy 
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^  _1 _ o£2  f  ^ 
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and,  after  the  appropriate  substitutions  for  ,  E„ ,  and 


d^G\ 


have  been  made. 
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1/k 

Ma  a  ^  1 

%P  _5/^  ^ 

X 
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The  surface  friction  is  thus  independent  of  the  viscosity  of  the  fluid 
and  varies  inversely  as  the  5/4  pover  of  the  distajice  from  the  slit. 

The  shearing  stress  is  independent  of  the  viscosity  of  the  fluid 
not  only  at  the  surface  but  also  between  any  two  adjacent  layers  in  the 
flow.  For  the  flow  from  a  slit  without  the  presence  of  the  plate  (refer¬ 
ence  8),  th?  shearing  stress  between  any  two  adjacent  layers  of  fluid  is 
also  noted  to  be  independent  of  the  viscosity  of  the  fluid.  The  shearing 
stresses  between  adjacent  layers  are  independent  of  the  viscosity  in  both 
the  problem  presented  herein  and  in  the  flow  from  the  slit  without  the 

plate,  because  in  both  cases  the  velocity  derivative  —  is  Inversely 

Sy 

proportional  to  the  viscosity  n  of  the  fluid. 

The  ratio  of  the  rate  of  flow  of  momentum  parallel  to  the  plate  at 
distance  x  from  the  slit  to  the  rate  of  momentum  flow  at  the  reference 
distance  a  from  the  slit  can  be  obtained  from  equation  (10)  as 

^  .  M  ^ 

The  rate  at  which  momentum  parallel  to  the  plate  passes  through  a  plane 
normal  to  the  plate  therefore  varies  inversely  as  the  1/4  power  of  the 
distance  from  the  slit  to  the  normal  plane. 

The  expression  for  |  can  be  obtained  as  follows: 


_  a,Ti  _ 
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The  width  of  the  disturbed  region,  given  by  y  for  f  =  £  .  thus 

^edge^ 

Increases  as  the  3/lj-  power  of  the  distance  from  the  slit. 


•  •  • 
•  •• 


Validity  of  Solution 

A  restriction  on  the  region  in  which  the  solution  is  valid  can  he 
obtained  hy  noting  that  for  the  houndhry— layer  eq^uations  to  he  valid  the 

value  of  ^  must  he  small.  By  use  of  equations  (26)  and  (28)  the  ratio 
becomes 


9I/3 
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J-/2  d| 
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3 _ /PV 

,V3^  V^a 


1/1+  1/12 

X  a 


Ma 

Large  values  of  — —  give  small  values  of  and  so  the  region  in 

PV^  u 

which  the  solution  Is  valid  is  enlarged  when  —  is  increased 

PV^ 

By  use  of  the  complete  equations  of  motion  instead  of  the  boundary- 
layer  equations  the  difference  between  the  static  pressure  on  the  plate 
and  the  static  pressure  very  far  above  the  plate  is  given  approximately 
at  a  given  value  of  x  by 


10  1  /  pv 


TT  2  _l/6  _l/2  I M, 

P'^max  j-  a  x  ^  ^ 
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and  the  pressure  gradient  in  the  direction  of  x  is  given  approximately 
by 


to*)/ 


1  ( PV^ 

_3/2\ 

X 


2/3 


When  it  is  noted  that  v  is  a  small  quantity,  the  approximation 
that  the  static  pressure  is  constant  throughout  the  entire  field  which 
was  made  in  equations  (l)  and  (3)  is  Justified  hy  the  solution.  All  the 
terms  neglected  in  the  houndary— layer  approximations  can  likewise  he 
shown  to  he  truly  negligahle  on  the  basis  of  the  obtained  solution.  All 


the  approximations  become  more  accurate  as  x 
increases. 


increases  and  as 


% 


because 


The  constant  a  seems  to  appear  in  a  rather 
Ma  /ai\^/^ 


where  a^  is  an  arbitrary 


is  the  momentum  corresponding  to  the  distance 
be  shown  that  a  change  in  a  does  not  change  u, 
vided  that  x  and  y  are  kept  fixed. 


arbitrary  manner  but, 

value  of  a  and 

from  the  slit,  it  can 
I,  or  Tq  pro- 


When  the  plate  is  removed,  the  flow  becomes  the  same  as  the  flow 
in  a  two-dimensional  laminar  Jet.  This  problem  has  been  treated  in 
references  7  and  8.  The  removal  of  the  plate  in  the  present  problem 


makes 


Equation  (20)  then  becomes 


d^G  „  d^G 
-  +  G  -  + 


d|- 


G  =  0 


0 


When  the  conditions  that 


at  ^  =  0  ajid  that 


"loo 

Jo 


are  used,  the  solution,  vhlch  in  this  case  is  obtained  in  closed  form, 
is 


Gr  =  tanh  - 
2 


(3^^) 


Because  the  momentum  Mq^  is  the  momentum  for  one  side  of  the  plate  in 
contrast  to  the  momentum  used  by  Bickley  Mg  (reference  8)  vhich  was 
the  total  momentum  flow  from  the  slit  and  therefore  eq^ual  to  twice  Mq^, 
the  result  |  =  21^  is  obtained.  When  this  relation  between  the  value 
of  ^  of  the  present  analysis  and  the  value  of  |g  in  reference  8  is 

used  together  with  equation  (3^),  exactly  the  same  relations  are  obtained 
for  all  the  quantities  of  Interest  as  were  obtained  by  Bickley 
(reference  8). 


CONCLUSIONS 


The  laminar  flow  of  a  slightly  viscous  incompressible  fluid  that 
issues  from  a  slit  and  passes  over  a  flat  plate  is  investigated  in  a 
region  far  enough  from  the  slit  for  the  boundary— layer  equations  to  be 
valid.  The  results  obtained  are  that  the  velocities  parallel  to  the 
plate  vary  inversely  as  the  1/2  power  of  the  distance  from  the  slit, 
that  the  velocities  perpendicular  to  the  plate  vary  inversely  as  the 
3/^  power  of  the  distance  from  the  slit,  that  the  width  of  the  disturbed 
region  Increases  as  the  3/^  power  of  the  distance  from  the  slit,  that 
the  surface  friction  is  independent  of  the  viscosity  of  the  fluid  and 
varies  inversely  as  the  5/^  power  of  the  distance  from  the  slit,  that 
the  rate  at  which  momentum  parallel  to  the  plate  passes  through  a  plane 
normal  to  the  plate  varies  Inversely  as  the  l/L  power  of  the  distance 
from  the  slit  to  the  normal  plane,  and  that  the  quantity  of  fluid 


•  ••• 
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passing  through  a  plan©  normal  to  the  plate  per  unit  time  varies  directly 
as  the  l/h  power  of  the  distance  from  the  slit  to  the  normal  plane. 


•  •  • 
•  •  • 
•  •• 
•  • 
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Figure  1.-  Flow  from  a  slit  at  0  into  the  XY -plane  with  a  plate  along  OX. 
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